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Abstract. — Wc obtain a classification of stationary, 7+-regular, non-degenerate 
and analytic electro-vacuum space-times in terms of Weinstein solutions. In par- 
ticular, for connected horizons, we prove uniqueness of the Kerr-Newman black 
holes. 
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1. Introduction 

We address the following celebrated and long-standing conjecture: 

Conjecture 1.1. — Let (^#, g, F) be a stationary, asymptotically flat, electro- 
vacuum, four- dimensional regular space-time. Then the domain of outer commu- 
nications ((^#oxt)) is either isometric to the domain of outer communications of a 
Kerr-Newman space-time or to the domain of outer communications of a (standard) 
Majumdar-Papapetrou space-time. 
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Arguments to this effect have been given in the literature [5,21,26] (see also [22, 
36]), with the hypotheses needed not always spelled out, and with some notable 
technical gaps. The aim of this work is to give continuation to the project initiated 
in [13], where the vacuum case was considered, and obtain a precise classifica- 
tion of such electro- vacuum solutions in the class of analytic space-times with non- 
degenerate event horizons, providing detailed filling of the gaps alluded to above. 

As usual, in mathematical Relativity, part of the challenge posed by a conjecture 
is to obtain a precise formulation. In the case of the "no-hair" conjectures this 
non-linearity lies in the notion of regularity and it is our opinion that the non- 
existence of a precise formalization for this concept has led to the enclosure of 
"hairy" assumptions and technical difficulties which has made the state of the art 



concerning this problem difficult to asscss|^j So, we start exactly by collecting 
our technical assumptions in the following (we refer to Section [5] for the necessary 
intermediary definitions) : 

Definition 1.2. — Let (^%,q) be a space-time containing an asymptotically flat 
end J^ext, and let K be a stationary Killing vector field on We will say that 
(^#, g, K) is I + regular if K is complete, if the domain of outer communications 
((^ext)) is globally hyperbolic, and if ((^# ox t)) contains a spacelike, connected, 
acausal hyper surf ace 5? D J^cxt? the closure 5? of which is a topological mani- 
fold with boundary, consisting of the union of a compact set and of a finite number 
of asymptotic ends, such that the boundary dy := 5? 7 \y is a topological manifold 
satisfying 

(i.i) oy c s+ ■.= d«^f C xt)) n i+M'cxt) , 

with dS" meeting every generator of S + precisely once. (See Fiaure [Ll\ ) 

Needless to say, all these conditions are satisfied by the Kcrr-Ncwman and the 
Majumdar-Papapetrou solutions and, in particular, by Minkowski and Reissncr- 
Nordstrom. For a detailed discussion of the previous definition and an alternative 
formulation of Conjecture 11.11 with regular replaced by a specific set of weaker con- 
ditions see [13]. 




Figure 1.1. The hypersurface y from the definition of 7 + -regularity. 



WAn illustrative example is given by the product structure 113. 151 1 that, although clear for 
Minkowski with the usual R X C/(l) action by isometries, seems far from obvious in the gener- 
ality required. Other examples are the regularity of the horizon and the asymptotic behavior of 
the relevant harmonic maps. 
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In this work we establish the following special case of Coniccture ll.il 

Theorem 1.3. — Let (^,Q,F) be a stationary, asymptotically flat, I + -regular, 
electro-vacuum, four- dimensional analytic space-time, satisfying (|2.7p and (12.111) . 
// each component of the event horizon is mean non-degenerate, then ((^# xt)) is 
isometric to the domain of outer communications of one of the Weinstein solu- 
tions of Section [U In particular, if the event horizon is connected and mean non- 
degenerate, then ((^ ex t)) is isometric to the domain of outer communications of a 
Kerr-Newman space-time. 

It should be emphasized that the hypotheses of analyticity and non-degeneracy 
are highly unsatisfactory, and one believes that they are not needed for the con- 
clusion. Note that by not allowing the existence of the "technically awkward" [6] 
degenerate horizons we eliminate extreme Kerr-Newman as well as the Majumdar- 
Papapetrou solutions from our classification. One also believes, in accordance with 
the statement of Coniecture ll.il that all solutions with non-connected event hori- 
zon are in the Majumdar-Papapetrou family; consequently one expects all other 
(non-connected) Weinstein solutions, and in particular the ones referred to in the 
previous result, to be singular. We postpone further discussion of this issues to the 
final section. 

A critical remark comparing our work with the existing literature is in order; 
we focus on those points that do not generalize immediately when passing from 
pure to electro- vacuum. First of all, the famous reduction of the Einstein-Maxwell 
source free equations to a singular harmonic map problem requires the use of Wcyl 
coordinates. The local existence of such coordinates has been well known for some 
time now, but global existence has, to our knowledge, cither been part of the ansatz, 
usually implicitly, or based on incorrect or incomplete analysis. The main reasons 
for this unsatisfactory situation resides in the existing proofs of non-negativity of 
the area function (|3. 10|) in ((^ C xt)), and existence of a global cross-section for the 
Ix U(l) action again in ((^# xt))- In [13] we solved this problems for vacuum 
and in Section [3] we present the necessary adjustments to extend this global exis- 
tence result to the electro- vacuum scenario. Also, no previous work known to us 
establishes the asymptotic behavior, as needed for the proof of uniqueness, of the 
relevant harmonic maps. More specifically: the necessity to control the behavior at 
points where the horizon meets the rotation axis, prior to [13], seems to have been 
neglected; at infinity, which requires special attention in the electro-vacuum case, 
part of the necessary estimates were imposed as extra conditions, beyond asymp- 



totic flatness; also, an apparent disregard for the singular character, at the axis, 
of the hyperbolic distance (|4.14p between the maps, even at large distance, appears 
to be the norm. A detailed asymptotic analysis is carried out in Section [5] 

We also note that a considerable part of the foundations of the theory underlying 
the desired classification depend exclusively on stationarity, / + -rcgularity and the 
null energy condition. Again, this work was carried out in [13], where various results 
were established under conditions weaker than previously cited, or were generalized 
to higher dimensions; this is of potential interest for further work on the subject. 



( 2 >See, for example, Theorem 2 in [27]. 
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2. Preliminaries 

An electro-vacuum space-time is a triple (^#, Q, F), assembled by a (n + 1)— 
dimensional Lorcntzian manifold (^#, g) endowed with a 2-form i* 1 , that satisfies 
the source free Einstein-Maxwell field equations 

(2.1) Ric-iR = 2T F , 

(2.2) F = dA , 

(2.3) d*F = 0, 

where Ric is the Ricci curvature tensor of the metric g, R its scalar curvature and 
Tp is the energy-momentum tensor of the electromagnetic 2-form F, 

(2.4) T F (u,v) :=$(i u F,i v F)- -\F\ 2 q{u,v) . 

A space-time will be said to possess an asymptotically flat end if ^ contains 
a spacelike hypersurface o5^ ex t diffcomorphic to R™ \ B(R), where B(R) is an open 
coordinate ball of radius R, with the following properties: there exists a constant 
a > such that, in local coordinates on c5^ e xt obtained from R" \ B{R), the metric 
7 induced by g on J^ext, and the extrinsic curvature tensor Kij of S^ cx t, satisfy the 
fall-off conditions 

(2.5) 7<i - % = O fc (r-«) , ^ = Ok-^r- 1 -*) , 
for some k > 1, where we write / = Ofc(r Q ) if / satisfies 

(2.6) d kl ...d k J = 0(r a - e ), 0<£<k. 

In connection with the field equations we also request the following decay rate for 
the electromagnetic potential 

(2.7) A M = O fc (r-°). 

A Killing vector K is said to be complete if for every p £ J% the orbit </>t [K] (p) 
of K is defined for all t £ R, i.e., if (the flow of) if generates an action of R by 
isometries; in an asymptotically flat context, K is called stationary if it is timelike at 
large distances. The exterior region -^> x t and the domain of outer communications 
^ext)) are then defined as 



(2.8) = / + ( u t &(^cxt) ) n /-(u t <M^cxt)) , 

with the event horizon being 

(2.9) #:=S((^e Xt )) ; ^ := /^ext) n <T . 

One expects stationary electro-vacuum space-times to be static or stationary- 
axisymmetric: static meaning that the stationary Killing vector is hypersurface- 
orthogonal, i.e., 

(2.10) dK b A K b = , 

where K b = g(K, •), and stationary- axisymmetric corresponding to the existence of 
a second complete Killing vector K^, which together with the stationary Killing 
vector if( ) := K generate an R x U(l) action by isometries. In connection with 
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the field equations we require the electromagnetic field to be invariant under the 
flow of the relevant Killing vectors 

(2.H) &k w F = 0. 

In the stationary and asymptotic flat scenario one is able to choose adapted 
coordinates so that the metric can, in a neighborhood of infinity, be written as 

(2.12) = -V 2 (dt. + Mf) 2 + j ij dx i dx j , 

—0 — T 

with 

(2.13) K = d t =^ d t V = d t 6i = dtjij = ; 

since we are also assuming electro- vacuum we get the following improvement of the 
original decay rates [7, Section 1.3], 

(2.14) ll3 - 5 l3 = OUr- 1 ) , 9i = O^ 1 ) , V-l = O^ 1 ) , 
and 

(2.15) A li = oo (r~ 1 ) , 

where the infinity symbol means that (|2.6p holds for arbitrary k. 



3. Weyl coordinates 

On a region charted by Weyl coordinates the source free Einstein-Maxwell equa- 
tions simplify considerably. It has been for long expected and recently showed in [13] 
that such global chart is available away from the axis of a stationary and axisym- 
metric vacuum domain of outer communications. In fact the role of the vacuum 
field equations in the referred analysis they imply the orthogonal integrability con- 
ditions (|3.ip and allow us to show that, whenever defined, the squared root of the 
area function (|3.10|) is harmonic with respect to the orbit space metric is fulfilled 
by the electro- vacuum field equations. 

The first of these well known results, which neither requires K(o) to be stationary, 
nor -fQi) to be a generator of axisymmetry, generalizes to higher dimensions as 
follows (compare [4]): 

Proposition 3.1. — Let(^,Q,F) be an (n+1) -dimensional electro-vacuum space- 
time, possibly with a cosmological constant, with n — 1 commuting Killing vector 
fields satisfying 

&K M F = , n = 0,...,n-2. 
If n — 2 of the zero sets stf^ := {p G \ Kr^\ p = 0} are non-empty thenv- 3 ' 

(3.1) dK (fl) A K (0) A . . . A K {n _ 2) = , V/x = 0, . . . , n - 2 . 

Proof. — To fix conventions, we use a Hodge star defined through the formula 

aA/? = ±(*a,/3)Vol, 

where the plus sign is taken in the Riemannian case, minus in our Lorentzian one, 
while Vol is the volume form. The following (well known) identities are useful [22]; 

(3.2) **9 = {-l) s(n+1 - s) - 1 , V6»gA s , 



( 3 'By an abuse of notation, we use the same symbols for vector fields and for the associated 
1-forms. 
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(3.3) i x *0 = *{6 AX) , V6»eA s , X e A 1 . 

Further, for any Killing vector AT, 

(3.4) [Sf K ,*] = 0. 

The Leibniz rule for the divergence 6 := *d* reads, for 9 E A s , 

S{9 A K) = *d * (6 A K)^ * d(i K *6) = *(3? K *9-i K d*9) 

113 10 * *J? K 6 - « A .(_i)(»+i-+i)(«+i-(«+i-+D)-i * * d * 

= (-lJ'^-'J-^jffl - (_l)»(n+l-»)-n+l * ^ A 

= {-i) a{ - n+l -^- l ^ K e + (-i) n+1 se a /v . 

Applying this to 9 = dK one obtains 

*d * {dK A A") = -& K dK + (-l) n+1 <taif A AT 

= (-l) n+1 5dK A K . 

As any Killing vector is divergence free, we see that 

SdK = (-l) n AK = (-l)"2trV 2 A' = (-l) n+1 2i K Ric , 

where A is the Laplace-Beltrami operator. The assumed field equations (with 
mological constant A) imply 

2 

Ric= 2 Tp + Ag , 

n — 1 

from which 

*d * (dK A K) = (-l) n+1 (-l) n+1 2ijf(2T F + — - — Ag) A A' 

n — 1 

= 2 ^2 i K T F AK + ^-j-A A'j A AT = 4 t A -T F A AT . 

Letting a := iif-F for any vector field X we have 

a ■ ixF = - * (a A *i x F) = -(-1)" * (*ixF A a) 

= (-l) n+ H a * *F = (-l) n+1 (-l) n+ H a i x F 
= -F(a,X), 



which inserted into (|2.4[) gives 

1 



i K T F = -i a F-^\F\ 2 K , 



and consequently 



*d*(dK AK) = —A(i a F + —\F\ 2 K) A K = -4i a F A K = 4 K A i a F 
Meanwhile, since (modulo sign) 

i a K = ± * (A" A *a) — ± * (K A *%kF) = ± * (K A *in * *F) 
= ± * [K A * * (*F A K)) = ± * [K A *F A K) = , 
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for (3 :— Ik * F £ A" 2 , we have 

*(a A 0) = (-i)ix(«-2) + A q) = (_i)n-2 ia * ^ 

= (-l)"- 2 i Q *i K *F= (-l) n - 2 i a * *(F A K) 
= (-l) n - 2 i a (~l) 3( - n+1 ^~ 1 F A /v 
= -(i a F A A' + (-1) 2 A A i a K) = -{i a F A A' + 0) 
= K A i Q A 

which leads to the significant 

(3.5) d* (dK A K) = 4a A/3 = ii K FAi K *F . 

Now, for any two commuting Killing vectors and an arbitrary differential form we 
have 

[-Sfjr j>) ,ijr Ci/) ]0 = ^K (li) (iK (u) d) - iK M (-&K M 0) 

= se Kw [6{K {U) ,...)}- {JZ Kw 6){K {v) , . . .) 

= (J? K( J)(K {U) ,...) + 0(JSfjc w if W ,-..)- (^k { J){K [v) ,...) = , 
giving us the commutation relation 

(3-6) [# w ,A- (l/) ] = [JSfjc M ,**jf M ] = , 

from which it follows that 

dF(K {u) ,K M ) = di Kw a [v) = -i K{ii) da {v) +Sf K<Jl) a^) 

= -%k w (-*K (v) dF + Sg Kw F) + i Kw Sf Kw F = 0, 

where we used the fact that F is exact and invariant under the flow of this Killing 
vectors. By the hypothesis on the zero sets, for any pair /i ^ is, we may take 
■^(p.) 7^ We then have F(K^, K^)\^, . = and consequently 

(3.7) F(K w ,K (u) )=0 , V/i,i/G {0,...,n-2} . 
A similar computation leads to 

(3.8) **(,.) **oo * F = , V^i, ^ e {0, ...,n - 2} . 
Now, let ijJ(fi) be the /x'th twist form, 

u ifi) := *(dK w AJf M ) . 

The identity 

^k w ^(v) = -&K W * (dK M A AT (ly) ) 

= *(&K M dK (v) + dK (v) A &K M K iv) ) = , 

together with 

-^Oi) • ■ ' ^(^) ^Wi)) = ^Oia) ■ ■ ■ ^^-iJ-^O^^Wi) = ' 

and Cartan's formula for the Lie derivative, gives 

(3.9) d{i K(n) . . . i^ )W(w+l) ) = ■ ■ ■ ^.odwcw+i) • 
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Wc thus have 

d*{dK (jt0 ) A K (fM]) A ... A A( Mn _ 2 )) = d(i K( ^_ 2) . -.i K{n) * (dff ((tio) A A (w)) )) 

-2) 



(-!)" 2 ^(,„_ 2 ) •■■^(, 1 ) (iw (Mo) 



1 13.51 , , „ . 

V ' (M2) 

= 4(— 1)%^, . , 

(F(^ (Mo) ,A' (Ail) )/3 (Mo) -a (Mo) Ai K(fll) iif ((i0) * F) 

So the function *(dFT( Mo ) A A( Mo ) AKr^ A . . . A-K> Mn _ 2 )) is constant, and, as before, 
the result follows from the hypothesis on the zero sets. 

□ 

Noting that a globally hyperbolic, stationary and asymptotically flat domain of 
outer communications satisfying the null energy condition is necessarily simply- 
connected [16, 19,20], in view of the previous result Theorem 5.6 of [13] translates 
to the electro-vacuum setting as: 

Theorem 3.2. — Let (^#,g,.F) be a four- dimensional, I + -regular, asymptotically 
flat, electro-vacuum space-time with stationary Killing vector K^ and periodic 
Killing vector Km, jointly generating an M X U(l) subgroup of the isometry group 
of(^£,Q). If ((^# xt)) is globally hyperbolic, then the area function 

(3.10) W:= -det ( S (K M ,K( V) )\ 

\ / fi,u— 0,1 

is non-negative on ((^# cx t)), vanishing precisely on the union of its boundary with 
the (non-empty) set {Q(Km, Km) = 0}. 

Away from points where A Km vanishes, which according to [13, Corollary 



3.8] correspond, in a chronological 4 ^ ((^# xt))j exactly to axis points 

(3.11) ^:={?€^ I Km\ q = 0}, 

there is a well defined and differentiable local cross-section for the Rx U(l) action. 
We can endow this cross-section with the orbit space metric 

(3.12) q(Z 1 ,Z 2 )= S (Z 1 ,Z 2 )-h» v &(Z 1 ,K (ll) )0(Z 2 ,K ( „ ) ) , 

whenever h^ v := q(K^,K^) is non-singular. The established orthogonality con- 
ditions allow us to identify, at least locally, the previous orbit space structure with 
a 2-surface orthogonal to the Killing vectors, provided by (|3.ip . endowed with the 
induced metric. From this and Theorem (|3.2[) we see that q is well defined and 
Riemannian throughout ((^# C xt)) \ it is then well known [36] that 



(3.13) A q VW = 0, 

whenever W is non-negative and q is Riemannian, which again is the case within 

((^ext)) \ Srf- 



( 4 'No closed timelike curves allowed. 
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According to the Structure Theorem [13], / + -rcgularity allows for the decompo- 
sition 

(3.14) «^ext» (I / + (.-#ext) =lxy, 

with Kn) tangent to 5? ', a simply-connected spacelike hypersurface with boundary 
which is an asymptotically flat global cross-section for the action generated by the 
stationary vector. We are now allowed to use the classification of circle actions 
on simply-connected 3-manifolds of Orlik and Raymond [29,31] to obtain a global 
cross-section for the K x U(l) action in ((-# e xt)) \ s& > Then, by (|3.13j) and relying 
on the results of [10], while disallowing the existence of degenerate horizons, we are 
able to undertake an analysis leading to 

(3.15) ((.-# cx t)) \^«lxS 1 xI+xl, 

while showing that this diffcomorphism defines a global coordinate system (t, (p, p, z) 
with 

(3.16) K {0) = d t , K {1) =d v and p = Vw . 

After invoking ()3.1j) once more, the desired global expression for the space-time 
metric in terms of Weyl coordinates 

(3.17) = -p 2 e 2X dt 2 + e- 2X (dip ~ wdtf + e 2u {dp 2 + dz 2 ) , 
follows, with 

(3.18) u = O fe _ 4 (r-- 1 ) , r = ^p 2 +z 2 -> oo . 



4. Reduction to a harmonic map problem 

The electro- vacuum field equations (|2.ip - (|2.3p and simple-connectedness of ((^# C xt)) 
guarantee the global existence of the following potentials: 

(4.1) dx = iK {1) F , dip = iKm * F and dv = u — 2(xdip — ipdx) , 
where 

(4.2) W :=*(dtf{ 1) Air{ 1) ), 

is the axial twist form. As discussed in detail in [36], when a global representation 
in terms of Weyl coordinates like (|3.17[) is allowed, the space-time metric is uniquely 
determined by an axisymmctric harmonic map 

(4.3) $ = (A, V ,x,^):R 3 W^H2 , 

here = {(0,0, z) z g I] and H^. is the 'upper half-space model' of the 2- 
dimcnsional complex hyperbolic space, i.e., K 4 with metric given by 

(4.4) ds 2 = d\ 2 + e 4X (dv + X dip - Hx f + e 2X {d X 2 + dip 2 ) . 

The metric coefficient A is part of the harmonic map and the remaining unknowns 
of the metric can be determined from $ by considering the unique solution (w, u) 
of the set of equations 

(4.5) d p w = -e iX p uj z , d z w = e 4X p v p , 
(4.6) 

dpU—d p X = p 



(d p X) 2 - (d z Xf + -/ X {lo 2 p - lu 2 ) + e 2X ((d pX ) 2 - (d zX ) 2 + (^) 2 - (9^) 2 



(4.7) d z u -d z X = 2p 



d p X d z X + \e 4X uj p co z + e 2X (d pX d zX + d z ?P) 
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that go to zero at infinity, and where we write uj a := u>(d a ) for a 6 {p, z}. 

4.1. Distance function on the target manifold. — The criteria for uniqueness 
of harmonic maps used in this paper (see Theorem 16.11 and compare [12, Appendix 
C]), is stated in terms of the pointwise distance between the maps. For the 'disk 
model' of Hp the distance between two points z = (z\, z-i) and w = (w\, W2) is given 
by [35, eq 55, pg 26] 

(4.8) cosh(d) = 



l-\z\\l-\w\ 2 



To obtain the distance function for the 'upper half-space model' we will use the 
isometry between the two referred models presented in [35, Appendix]: first we 
perform the coordinate transformation 



1 — xi 2x 



z\ = — , z 2 



2 



1 + Xi 1 + X\ 

with analogous expressions for w% = 1/^(2/1,2/2) to obtain 

2 |xi +2/1-2^2 2/2 1 



|1 — Z±Wl — 22^2 1 



l + a?i [l + 2/i I 



then we take 



so that 



= + and = -H±M 



(4.9) cosh(d) = i |xi + 2/1 - 2«2 2/2I e Al+A2 ; 
and finally, by writing 

(4.10) xi = e~ 2Al + Xi + ^1 + 2 wi and x 2 = Xi + #1 , 

with similar expressions for 2/j = Di(M, ^2, X2, ^2), we see that the distance function 

satisfies |w 

(4.11) 

COsh 2 (d) = 1 e 2(A 1 +A 2)(e -2A 1 + g -2A 2 + ^ _ ^2 + ( ^ _ ^ )2) 2 



+ e 2(Al+A ^ 2 - Vl - Xli'2 + X21P1) 2 

= \ {e" Al+A2 + e Al - Aa + e Al+A2 (xi - X2) 2 + e Al+A2 (^i " ^) 2 } 

+ e 2 ^ +x ^ {(« 2 - «i) + (X2^i - Xi^ 2 )} 2 , 
or in an apparently more intrinsic way 
(4.12) 



2 



\ 2 

2 1 (,h J, \2 I 



cosh 2 (d) - 1 ' l 02 ^' 9 ^ ■ I3i(9<p,d v ) , (xi - X2) + (V>i - ^2) 



4 [ y 0i(9 v ,5 v ) y 82(6?, d v ) ^31(^,^)^02(9^,^) 
f («2 - vi) + (x2-0i - XlV^) 1 



( 5 ) By taking Xi = ^1=0 wc see that this distance function is related to the one used in the 
vacuum case [13, Section 6.5.1] by d = 2df,. This discrepancy has its genesis in an analogous 
relation between the line elements of the different disk models used. 
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Figure 5.1. The quotient space and its double. 



It will also be helpful to use the usual rescaling 

(4.13) Ui = Xi+hip, so that Qi(d v ,d v ) =p 2 e~ 2U > = e - 2X > , 

from which we get our final expression for the distance in the 'upper half-space': 
(4.14) 

cosh 2 (d) = 1 {e u *- u > + e~ u ^ u * + p-*e u * +u * ( X i - X2) 2 + p~ 2 e u ^ (Vi - ^ 2 ) 2 } 2 
+ {p- 2 e u >+ u *(v 2 - Vl ) - p- 2 e u ^( X i^2 - X2^i)} 2 ■ 



5. Boundary conditions 

5.1. The Axis. — From now on we will be controlling the distance, as given by 
any of the formulae in the previous section, between the harmonic maps arising from 
two / + -regular, stationary-axisymmetric and electro- vacuum space-times 
i = 1,2. We will start by showing that 



(5.1) d($i, $ 2 ) is bounded near stf n ((^# e xt)) ■ 

In this section we will be working with the following coordinate systems: isothermal 
coordinates (xi, z{) globally defined in the doubling across the axis of the orbit space 
of an appropriate extension of the U(l) action to the manifold obtained by the 



addition of 3-discs to every connected component of dS^i "canonical coordinates" 
(p, z) of the half plane Rg~ x R, which is the image of each (physical) orbit space by 
the map defined by {x l ,z i ) i-> (p l (x l , Zi), z l (x l , zi)). 

Let <j> s be the flow generated by the axial Killing vector Kny In the doubling of 
the orbit space the isothermal coordinates satisfy 

x o X = —x and z o X = z . 

Then, invariance of a function (x, z) 1— > /(x, z) under the axial flow, which is the 
case for the fields v,x and ip, implies that the function x t— > f(x, z) is even for all 



( 6 )The resulting space is diffeomorphic to M 2 , see Figure [5TT1 and for more details concerning this 
construction see [13, Section 6]; also, the fact that 7+ -regularity, stationarity and the null energy 
condition imply spherical topology for the connected components of the cross-section of the event 
horizon QS^i follows from [16]. 
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z. In this case, if / is C 2 , Taylor expanding on x, from the axis, gives 

(5.2) f(x, z) = /(0, z) + l^l( c (x),z)x 2 , |c(f )| < \x\ • 

Now fix a point in si n ((^# C xt)) and by rescaling I assume it lies at the origin. Sup- 
pose also that / = /o := /(0, 0) along the connected component of ^ H ((^Cxt)), 
in ((^# C xt)), containing (0, 0); this is clearly the case for all the functions appearing 
in 14.111 and it also implies that we can realize the aforementioned extension of the 
doubling of the orbit space to M. 2 while preserving the constancy of / along the 
extended axis near the poles, i.e., near the points where the axis meets the event 
horizon. Then (|5.2|) implies 

(5.3) \f(x,z)-f \ <Cx 2 near (0,0) e si n ((.^ cx t)) . 

We will need bigger control over the functions e~ 2X — To this end let {x, y, z} 
be Gaussian coordinates along the axis, in the extension of J? 7 , with si = {x — y = 
0} and for which Km = xd y — yd x (see [10, pg 5] and compare with (|5.23| ). For 
any path with initial velocity transverse to si we have 

(5.4) S7^(o)K {1) \ x=y=0 = V 7 < 9i (xe?„ - yd x )\ x=y=0 = ^ x d y - -y y d x , 

and consequently g(\7^ {0) K {lh \7^ {0) K {1) ) = {l x ) 2 + {"i v ) 2 ^ 0. Since V^0(K^, Km) = 
2g(V \j_Km, we see that the gradient of g vv vanishes at the axis and 

V M V„ fl^d), -Sr ( i))U = 2 (V M A' (1) ,V^ (1) )U ■ 
Taylor expanding along 7 yields 

B v <p 7(s) = (fl(V^(o)-RT(i), V m K w ) +0(s))s 2 , 

from which it follows that for any path transverse to si and small s 

(5.5) C~ l s 2 < W 07(a) < Cs 2 . 



We will need to consider two separate cases. First, fix, in each space-time, a 
point belonging to si n ((-# e xt)) and rescale all the previous coordinate systems 
so that each of the fixed points corresponds to its respective origin and ^^(0,0) = 
(0, 0). At these points, since there the boundary of the orbit space is analytic, the 
function pi = p^(xi,Zi) may be extended analytically across the origin, therefore, 
as an immediate consequence of (|5.5p we get control over the first terms appearing 



111 



6H 



(5.6) e v '- v * = W^HH ^ \ ^ c near (°>°) e ((^ ext )) • 

Since the Xi's and the tp^s are all bounded near the origin our goal gets reduced to 
showing that 

(5.7) p- 2 e Ul+U2 (h - f 2 ) = 0(1) near (0,0) e sin ((^ cxt )> , 

when /1 = Xi) V , i)' y l)XiV'2 and f 2 = X2,ip2,V2,Xii>\, where by this we mean that 
if, for example, we set f\ = xi then f 2 = X2- 

Let us start with fi = xi, V"i) v i and /2 = X2i ^27 ^2- Each /, is invariant under 
the respective axial flow and constant along each connected component of sii n 
((^#ext)}) so if wc impose /i(0, 0) = /2(0, 0) = /o, which is always achievable if the 
space-times Qi) have the same set of masses, angular momenta and charges (see 
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Section [6] and [36, Section 2.3]), we sec that (|5.3|) holds and using (|5.5[) and (|5.6[) 

we get 



(5.8) 



-2„(7i + !7 2 



(/i-/a)| = 



/1-/2 



< 



l/i -/o| + |/ 2 -/o 



l/i - /o| Bi(9 9 ,d v ) |/ s 



-/o| Q2(d< p ,d lp ) 



< 



< 



C~ x x\ 



32(d v ,d ip ) 
:C 6 



C4X2 



1a2 



C near (0, 0) e ^ H ((^ ext )) 



We take the chance to stress the fact that the previous argument does not apply 
to the fields xii/^ and XaV'i since these products involve functions originating from 
different space-times and therefore only make sense as functions of (p, z) for which 



{<) 



estimates like 15.31 are not available a priori 
To bypass this problem we write 

XlV>2 - X2^1 = (Xl + X2)(V^2 - Ipl) + XlV'l - X2^2 • 

Since Xi + Xi is bounded, to control the first term we just need to take fi = ipi as 
before. Setting fx = xiV'i an d fi = Xi^i we see that the previous argument still 
applies as these are also axially symmetric functions which are constant along the 
axis components. The desired result follows. 



To finish the proof of boundedness of (|4.14j) near the singular set stf n ((^#ext)) 
we still have to analyze what happens near points where the axis meets the horizon. 
Choose such a point in each space-time and, without loss of generality, assume that 
these are 'north poles' which, as before, lie at the origin of the coordinate systems 
(sbi, Zi), and satisfy (p, z) = ^(0, 0) = (0, 0). 

As already mentioned, a careful extension of the doubling of the orbits spaces 
validates (|5.3p in a neighborhood of these 'north poles', but, on the other hand, the 
Pi's are now non-differentiable at such points and (|5.6[) no longer holds. Nonetheless, 
if we are able to control e Ui ~ Uj by other means, then the inequalities established 
in (|5.8|) extend to the case under consideration and boundedness of the distance 
near the axis follows. This problem, which is in fact the major difficulty that arises 
in the analysis of the boundary conditions of this axisymmetric harmonic maps, has 
been recently overcomed for the vacuum case [13, Section 6.5.1] by obtaining the 
following uniform estimate 



(5.9) 



U = In 



\jz + yfz^Tp^ + O(l) near (0, 0) G si n S + 



from which the desired consequence immediately follows. This result, which requires 
this component of the horizon to be non-degenerate, extends to the electro- vacuum 
case immediately. 



("'in fact, extending pi and near this axis points by pi(— Xi, Zi) = —pii&i, Zi) and Zi(—Xi, Zi) = 
Zi(£i,Zi) shows that invariance under the axial flow implies that f(p,z) := f o ^/~ 1 (p,z) is an 
even function of p. Then, direct estimates in terms of p analogous to 1 15.31 1 may be obtained for all 
the fields and the presented procedure including 115. 81 may be bypassed. Unfortunately this is no 
longer possible near points where the axis meets the horizon as the pi are no longer diffcrentiablc. 
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5.2. Spatial infinity. — In this section we want to show that 

(5.10) lim d($ 1 ,$ 2 )=0, 

with d implicitly defined by (|4.14p . For this we will assume stationarity and asymp- 
totic flatness as given by the system of equations (|2.12[) - (|2.15|) . It turns out that 
the estimates provided by asymptotic flatness, even in the way just defined, seem 
insufficient to control the relevant fields; even in an adapted frame provided by 
the results of Section [3] integration of the defining equations (|4.1[) yields divergent 
logarithmic terms. Fortunately, in the stationary and electro-vacuum setting, the 
asymptotic analysis of Beig and Simon [3,32] provides relevant improvements of the 
initial decay rates by means of the expansion (15.151) . 

Let D and e^fc denote the covariant derivative and volume element of 7, the 
induced metric in .y ex t as in Section O A well known consequence of the source 
free Einstein-Maxwell equations ()2.ip - (|2.3|) and simple-connectedness of ((^ext)) is 
the global existence of functions r and a satisfying [24,32,36] 

(5.11) D t r = V'./ illAi. + 9 3 D k A Q ) , 
and 

(5.12) D t a = -V^e^DjOk + i(^d^ - *<%*) , 
where VP := A + it. If we introduce the Ernst potential 

(5.13) <*f = V 2 - + ia , 

and consider the complex valued fields £ and implicitly defined by 

(5.14) ¥ = , 



1 + ?? ' l + d 

then [32, eq 3.11] provides the following expansion for the vector S A := £) G C 2 
in terms of an arbitrary asymptoticaly flat coordinate system 

M A M A T k 

(5.15) * A = ^T + + Ooo(logr/r 3 ) . 

We note that the apparent discrepancy between the error term here with the one 
in the original paper comes from the fact that the result there is presented in 
adapted coordinates obtained from arbitrary asymptotically flat coordinates by a 
transformation of the form x % 1— > x % + O^logr). 
Using the identity -g^y = ^ - grgq^ we get 

y= l Z 4 = l+^-2tf =1 __2tf =1 _ 



1 + 1? l + i? l + i? V 1 + ?? 

^3 



= 1-2 tf-tf 



)2 



1 + 0, 

Inserting the ^-component of ()5.15j) into the last expression yields 

Sf = 1 - 2 ( — + + Ooo(logr/r 3 ) ) + 2 — + + 0O (l g r/r 3 ) ) + M (r 3 ) 

= 1 - 2^- - 2^- + 2^J- + 0oo(r -4) + Ooo( logr/r 3 ) . 

Noting that the topological restrictions imposed by asymptotic flatness imply that 
the imaginary part of M vanishes, = [32, Section IV], we write M = M 
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and by setting = Mk + iSt we get 

(5.16) V = 1 - 2^ + 2^ - 2^ - 2i^- + Ooo (logr/r3) . 
Consequently 

(5.17) (T = ^ = _ 2 ^ + 00 (logr/r 3 ). 
Similarly for "J we get 

C Ci? o? 2 

* = — — = C — - — = C - 6? + - — ■ 

Inserting (|5.15|) into the last expression yields 
* = C(l-^) + Ooc(r- 3 ) 

= I — + + Ooo (log r/r 3 ) ( 1- — ^3— + O00 (log r/r 3 ) J + Ooo(r 3 ) 

= — -2 — + —=5— + Ooo(logr/r 3 ) . 

As before 3M" = 0. So now, by setting = § and = Qfc + zA, we see that 



z r zr 



(5.19) r = 3* = ^ + 0oo (l g r/r 3) 



We have - = 2i(A)9iT - r<Mo) = O^fr 4 ) and using (12~T1 

and (|5.17[) we get 
(5.20) 

</ } l>/h = -V- A (D t a - - )) = A (2^- + Ooo (log r/r 3 ) 

With the exception of the already noted logr discrepancy in the error term, this 
is [3, eq 4.1, pg 1010] and so we get 

(5.21) 9i:=^ = 2e ijk ^- + O fc (log r/r 3 ) , 

Qtt r } 

where e [ljk] = e ijk with ei 23 = 1. 

5.2.1. The electromagnetic twist potential and the norm of the axial Killing vector. 

- Until now we have been working with a generic asymptotically flat coordinate 
system, but to estimate the electromagnetic twist potential v via the Ernst equa- 
tions (|4.5p and the results of the previous section wc will need to use adapted 
coordinates. So, let {t, p, p, z} be the Weyl coordinates as constructed in Section [3] 
and define the cylindrical type coordinates 



(5.22) 



x = p cos p 
y = p sirup 

A simple but noteworthy fact is that in this coordinate system we have 



(5.23) = d v = xd y - yd x . 

The estimates of the previous section will only be available to us in this coordinates 
if {t,x 1 } = {t, x,y, z} is an asymptotically flat coordinate system. This is in fact 
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the case. To see it note that in the orbit space {t — <p = 0} the identity [13, eq 6.9] 
yields 

(5.24) d p = {1 + Okif-^di + Oodf-^di (<p = 0), 

(5.25) d z =O oa {f- 1 )d i + {l + 00 (r- 1 ))d i (9? = 0). 

Recall that {x, z} are asymptotically flat isothermal coordinates (for the orbit space 
metric). Direct computations yield 2 xx \ v= q = 1 + O o( r ~ 1 )- a similar expression for 
Q yy and, using g pv = 0, also K2/ | v= o = 0. The defining decay rates are now obtained 
by flowing the previous estimates along the integral lines of the axial Killing vector. 
We illustrate this with an explicit calculation: 

$xy\(p=-<p = fKOApo)*^*, {4> Vo )*dy)\ V =Q 

= g(cos </5o<9x + shupody, — smtpodx + cos <fody) \<p=o 

= - sin ip cos ip q xx I v=0 (cos 2 ip - sin 2 ^ ) g xy | v=0 + sin ip cos ^ 9yy I v =o 

= 1+O 00 (r- 1 ) =0 =1+O 00 (r- 1 

= Ooc(r- 1 ) . 

So we have constructed asymptotically flat coordinates for which the following 
uniform estimate holds 

(5.26) 9 vv \ v =o = P 2 9yy\ v =0 = p 2 (l + Ooo^ 1 )) . 
As a nice consequence we get 

(5.27) e- 2U :=^- = l + O0 (r- 1 ), 

P 

from which we see that 

(5.28) e u * ±u > := (1 + O^r^Ml + O 00 (r- 1 )) ±1 = 1 + CUr" 1 ) - r ^ +00 1 , 
and our goal (|5.10p gets reduced to showing that 

(5.29) _Jmr ^ ~ ^ ? = Jim 

Vl - V 2 



lim 



2 



lim xii>2 - X2-01 = 



^/p 2 +z 2 -^+00 P 

It follows from (|5T2"Tj) and ([3"T7j) that 

3zt = 0^9 z = 0^S x = S y =0. 
So we set J := — S z and by using (|5.21[) with (|2.14p we get 

(5.30) 3vt \ v=0 = 2 J^ + Ooo (log r/r 3 ) , 
from which 

(5.31) v t| y =o = 6(pd y ,dt)\<p=o = 2J -^ + P°oo (log r/r 3 ) 
and therefore 

(5.32) iselj = ^ + I 0oo (io gr /r 3 ). 
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The Ernst equations (|4~5)) together with the estimates (|Q3)) , (|Q5)) , (jOT)) 
and (|5.48|) . that will be established in the next section, provide 



(5.33) 



d z v = — 6Jp 4 /r 5 + pOoa (log r/r 3 ) 
<9 p v = 6Jp 3 z/r 5 + p 2 Ooo (log r/r 4 ) 



Integrating this system by using the polar coordinates p = rcos8. z = r sin#. while 
imposing the standard condition 

(5.34) v{0,z) = 0, for z > , 
yields 

(5.35) « = 4J - I - ( 3p2 ~ z2 + 9] + P O x (log r/r 2 ) . 



2 r V r 

We note the following relevant relation with the total angular momentum as 
given by the Komar integral formula 

lim / *dKf-<\ = lim 2tt [ ijf„, * dK)-,^ 

+«= 16tt J {r=R} W ^ +00 16tt y {r=fi}n{¥ , =0} (1) 

1 



lim / *(^'(i) A -K'(i)) 

~ lim / (dv + 2 ( X # " V^x) ' 

8 R-» +cx) y{ r=fl }n{¥'=0} 



=0(r 



-\ lim («(0,i?)-«(0,-i?)) 



= J . 

We are now able to establish the electromagnetic twist potential part of (|5.29|) . 
For two twist potentials satisfying (|5.34|) we have 

(5.36) lim Vl ~ 2 V2 = . 

To take care of the asymptotic behavior of v near the axis we Taylor expand on 
p around a point (0, z), away from the poles, to get 

(5.37) v(p, z) = v(0,z) + d p v(c(p),z)p , \c(p)\ < p . 
Then, using (|5.33|) to obtain 

(5.38) d p v = p 2 00 (r- 3 ) , 
we conclude that for \z\ ^> and p < \z\ 

(5.39) v = v{0,z) + p 2 0(r- 3 ) . 

Finally, for two twist potentials that agree along the axis for both large positive 
and negative z we have, in the region p < \z\, 

(5.40) p- 2 (v 1 -v 2 )=0(r- 3 ) , 
and the desired result follows. 
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(5.42) 



5.2.2. The electromagnetic potentials. — Asymptotic flatness (|2.7[) together with (|5.18[) 
and (|5.19[) yield the desired improvement of the initial decay rates 

(5.41) = Ooo(r~ 3 ) . 

Now in the {x"} = {t, x, y, z} coordinates of the previous section we have 

dx ■= iK {1) F 

= F tlv dx»dx v (K {l) ,-) 

= F„ v (dx»{K {1) ) dx" - dx v {K {l) ) dx") = 2F^dx"{K (1) ) dx" 

= 2F^dx"{xd y - yd x ) dx" = 2F liV (x6$ - y6%) dx" 

= 2 {xF yu - yF xv ) dx" = 4 (xd[ y A^ - yd[ y A„]) dx" . 

With (|5.41[) we see that, in the orbit space {<p = 0} (where y = 0, x = p and 
dx = d p ), we have 

d P x\ v =o = 4pd [y A p] = pOooir' 3 ) 
d z x\ v =o = 4pd [y A z] = pO^r- 3 ) 
Imposing the boundary condition 

(5.43) x(0, z) = 0, for z > , 
integration yields 

(5.44) x| y =o = pOoo(r- 2 ) . 

Arguing as in the end of section 15. 2. II the equation in (|5.29p corresponding to the 
potentials Xi follows. 

To obtain a coordinate expression for dip it will be helpful to rearrange our 
preferred coordinate system and consider {x^} = {t, y, x, z}, then 

dip := i K(1) * F 

= ^F»"e^ Xa dx x dx°(K {1) , •) = F^e^ Xa dx x (xd y - yd x ) dx° 

= F^e^ixS* - yS x ) dx° = F^(xe„ vya -ye„ vxa ) dx° . 
Now, in the orbit space and away from the axis, we have (compare with (|3.17p ) 



/ 


Qtt 


P 1 Qt v 
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p~ 
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e 2u 


J 



(5.45) Q^\ v =o = 



therefore 

(5.46) det( ^| v=o ) = (stt^f ~ Q \\ e 4u = \{- P 2 )e iu = ~e iu 



so 

9 P V-| v =o - pe^ yx F^ = p(e tzyx F tz + e ztyx F zt ) = 2pe tzyx F tz 



2p^\dct( 9 ^)\F tz = 2pe 2u Q» t Q " z F llv 
2pe 2u Q^Q zz Fp Z = 2pe 2u e~ 2u ( 9 u F tz+& y t F yz ) 

{-^{d t A z -d z A t ) + 2^d {y A z ^ 



2 



P 

2^d z A t + 4 Stv d [y A z] = 2p(l + O^r-^Ao + pO^r" 5 ) 
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where in the last equality we used (|5.26[) . (|5.31[) and (|5.41[k also 
3 z ^| v =o = pe^yzF^ = 2pe txyz F tx 



-2p^/\det( 8 ^)\F tx = -2pe 2 VV^ 
^2pe 2u Q^Q xx F^ x = -2p( Q tt F tx + g^Fy*) 

-2p L^( 9t A x - 8 x A t ) + 2^d [y A x] } 
-2p(l + O oo (r- 1 ))9^ + pO oo (r- 5 ) . 



From (|5.18p we get 
(5.47) 



d p ^\ v=0 = -Qf#+pOoo(r- 3 ) 

Integrating as before while using a standard boundary condition provides 

(5.48) ^ = Q(-1+^) +pO oa (r- 2 ) . 

We note the following relevant and expected relation with the total electric charge 
given by the Komar integral 

(5.49) lim — - / *F = Q . 

+oo 47T J{ r=R } 

It should be now clear that (|5 . 29[1 follows. 

The results of this last two sections establish one of the significant missing ele- 
ments of all previous uniqueness claims for the Kerr-Newman metric: 

Proposition 5.1. — Let VE^ = (U i ,v i ,Xi,i'i), i = 1,2, be the Ernst potentials 
associated with two I + -regular, electro-vacuum, stationary, asymptotically flat ax- 
isymmetric metrics with non- degenerate event horizons. If Vi = «2, ipi = V'2 and 
Xi = X2 on the rotation axis, then the hyperbolic- space distance between and ^2 
is bounded, going to zero as r tends to infinity in the asymptotic region. 

6. Weinstein Solutions: existence and uniqueness 

In this section we construct axisymmetric Ernst maps 
$= (U,v, X ,il>) :K 3 W^Hjl, 

which are "close" to some reference maps, not necessarily harmonic, satisfying con- 
ditions modeled on the local behavior of the Kerr-Newman solutions. First recall 
the definitions of mass, angular momentum and electric charge of the fc-th black 
hole as given by the Komar integrals 

(6.1) m k := -— / *dK\ 0) , 



8 71 



s k 



(6.2) J„ : = -L J^K^ , 

(6.3) *°-=-hh; F - 

for some 2-sphere Sk whose interior intersects the event horizon exactly at its fc-th 
component. 
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We are now able to characterize the reference maps $ — (U, v, \, ip): 

1. The components f = v,X and ip are locally bounded, constant along each 
connected component of si \ S + = U^_ ^4 and we write f\^ k = These 
functions are normalized to satisfy /jv = 0. 

2. There exist Ndh > degenerate event horizons, which are represented by 
punctures ((p = 0, p = 0, z = bi), together with a mass parameter m, > 0. 
In a neighborhood of such puncture, containing only this component of the 
horizon, the map $ corresponds to the harmonic map of the (extreme) Kcrr- 
Ncwman solution parameterized by 

( \ i h+1 ~ ti \ 

(m l ,q l ) = (m,, ) . 

3. There exist N n ^ > non-degenerate horizons, which are represented by 
bounded open intervals (c^,cf) = Ii C si, with none of the previous b/s 
belonging to the union of the closures of the Ii. In a neighborhood of such in- 
terval, containing only this component of the horizon, the map $ corresponds 
to the harmonic map of the Kerr-Newman solution parameterized by 

(jJ,j,\j,qj) = (2 dz,v i+1 - Vi, . 

To retrieve the usual parametrization using mass, angular mommentum and 
charge one uses the known explicit formulas for Kerr-Newman (e.g., equations 
2.31. of [36]) together with the following relations [36, section 2.3.] 

(6.4) j, = h+k , 

(6.5) = //; + 2WjJj , 

where the auxiliary parameters are defined by Xj := f I dv,lj = J T xdip—ipdx, 
and w\i j =Wj, with w defined by (|3. 1 7[) . 

4. In a neighborhood of infinity the functions U, v, x and V* coincide with the 
components of the harmonic map associated with the Kerr-Newman solution 
with mass M '■= ^2 k m.k angular momentum J := Jk = wo/8 and electric 
charge Q := J2k 1 k = — V-'o/2, where the sums are taken over all the compo- 
nents of the event horizon. 

5. The functions U, v, x and ip are smooth across si \ (Ui{bi} Uj Ij). 

A collection {bi, rrii}f^, {Ij,v(cJ), vic^)}^^, and {ipk\k=o wm be called "electro- 
vacuum axis data". 

A map <!> satisfying condition 1.-5. above defines singular Dirichlet data [37, Def- 
inition 2] (compare [36, Section 2.4.]) with a target manifold with constant negative 
sectional curvature. We then have the following version of [37, Theorem 2] (com- 
pare [12, Appendix C] where the uniqueness claim is clarified, and [36] for a similar 
result stated purely in terms of axis data): 

Theorem 6.1. — For any set of electro-vacuum axis data there exists a unique 
harmonic map $ : R 3 \ si — ► whose distance, as given by (|4.14[) . from an ax- 
isymmetric map : R 3 \ si — > H^., not necessarily harmonic but with the properties 
1.-5. above, satisfies: 

(6.6) d($,$) G L°°(R 3 \0 , 
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and 

(6.7) d($, $) -> as r -> +00 . 

□ 

From an harmonic map $ : R 3 \ — > Hp one can construct a stationary and 
axisymmetric solution of the source free Einstein-Maxwell field equations [36, Sec- 
tion 4.1.]. Such (not necessarily I + -regular) space-times, arising from the harmonic 
maps of the previous theorem will be referred to as Weinstein solutions. 

7. Proof of Theorem [TT3l 

If <? + is empty we obtain Minkowski by an aplication of [7, Theorem 2.7]. Oth- 
erwise the proof splits into two cases, according to whether or not K is tangent to 
the generators of <f + . 

Rotating horizons: 

Suppose, first, that the Killing vector is not tangent to the generators of some 
connected component . Proposition 1.9 of [7] allows us to generalize [13, Proposi- 
tion 4.10] to electro- vacuum and then Theorem 4.11 together with the Remark 4.12 
of [13] show that the event horizon is analytic if the metric is; also, by (|5.4ip and 
Einstein's equations, G M „ = 2T fJtv = 0(r~ 5 ). So the Rigidity Theorem, as pre- 
sented in [8, Theorem 5.1], applies and establishes the existence of a K x U(l) 
subgroup of the isometry group of ', g). The analysis of Section 02 leading to the 
global representation (|3.17|) of the metric, is now available. As stressed through- 
out this paper, in this gauge, the field equations (|2.1[) - (|2.3p reduce to a harmonic 
map $ (|4.3|) . The analysis of the asymptotic behavior of such map, whose results 
are compiled in Proposition 15.11 shows that <& lies a finite distance from one of 
the harmonic maps associated to the Weinstein solutions of Theorem 16.11 and the 
uniqueness part of such theorem allows us to conclude; note that in the connected 
and non-degenerate setting the Weinstein solutions correspond to the non-extreme 
Kerr-Newman metrics. 

Non-rotating case: 

Now let us consider the case when the stationary Killing vector is tangent to 
the generators of every component of $ + . Following the procedure in [13, Section 
7.2], based on [30], we extend ((-# e xt)) to a space-time where each connected com- 
ponent of the event horizon is contained in a bifurcate horizon. Then, by [15] there 
exists an asymptotically flat Cauchy hypersurface for the domain of outer commu- 
nications, with boundary on the union of the bifurcate spheres, which is maximal. 
We are now able to conclude from Theorem 3.4 of [33] that ((.-# C xt)) is static. By 
taking in account the corrections presented in the proof of Theorem 1.4 of [13] wc 
can invoke, after relying on analyticity once more, the non-degenerate part of the 
conclusion of Theorem 1.3 of [9], yielding non-extreme Reissner-Nordstrom as the 
only non-rotating solution satisfying the remaining conditions of the desired result. 

8. Concluding remarks 

To obtain a satisfactory classification in four dimensions, the following issues 
remain to be addressed: 
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1. Analyticity. The previous versions of the uniqueness theorem required an- 
alyticity of both the metric and the horizon. As shown in the proof of Thco- 
rem ll.31 the latter follows from the former. This is a worthwhile improvement, 
as even C 1 - differentiability of the horizon is not clear a priori. But the hy- 
pothesis of analyticity of the metric remains to be removed. In this context 
one should keep in mind the Curzon solution, where analyticity of the metric 
fails precisely at the horizon. 

We further note that a new approach to Hawking's rigidity without analyt- 
icity [1,2] as yield significant breakthroughs in the vacuum case. According 
to A.D. Ionescu (private communication) the generalization of the results to 
electro-vacuum should follow by similar techniques. However, some problems 
still need to be settled, even for vacuum: the local claim requires a non- 
expanding horizon which we expect to be a consequence of 7 + -regularity and 
the results and techniques of [11, 13], but such claim requires checking; also, 
as it stands, the global result is restricted to near Kerr geometries. 

The hypothesis of analyticity is particularly annoying in the static context, 
being needed there only to exclude non-embedded Killing prehorizons [13, 
Section 5]. The nature of that problem seems to be rather different from 
Hawking's rigidity, with presumably a simpler solution, yet to be found. 

2. Degeneracy. The classification of black holes with degenerate components 
of the event horizon requires further investigations. We believe that the re- 
sults here go a long way to obtain a classification, in terms of Wcinstein 
solutions, of stationary, axisymmetric, rotating configurations allowing both 
degenerate and non-degenerate components of the horizon: the foundations 
are settled but we are still missing an equivalent of Proposition 15.11 Recall 
that in the static case a complete classification in terms of the Majumdar- 
Papapetrou and the Reissner-Norsdtrom families, with neither degeneracy or 
connectedness assumptions is already available by the work in [17] and ref- 
erences therein. In fact more is known in the degenerate class, since it was 
established in [14] that appropriately regular, J + -regular in particular, Israel- 
Wilson-Perjes Black holes belong to the Majumdar-Papapetrou family. 

It has been announced [23] that the question of uniqueness of degenerate 
black holes (with connected event horizon) has been settled. Unfortunately, 
that reference does not contain any new results, as compared to what had 
already been published in [13], or is contained in this work, and so, it is our 
belief that this problem remains open. Indeed, the existence of global Wcyl 
coordinates with controlled behavior at the singular set is assumed. In the 
non-degenerate case this issue was first settled for vacuum in [13], but the 
degenerate case appears to present serious technical difficulties, and requires 
further study. 

3. Multi Component Solutions. In agreement with the statement of Con- 
jecture II. 1[ one believes that all solutions with non-connected $ + are in the 
Majumdar-Papapetrou family. From what was said in the previous item, we 
see that it remains to show that non-static Weinstein solutions with non- 
connected horizons are singular; besides the already quoted result dealing 
with the Israel-Wilson-Perjes family, this has been established for slowly ro- 
tating black holes in vacuum by a regularity analysis of the relevant harmonic 
maps [25,34] and recent and promising results seem to have settled the prob- 
lem for two-body configurations, also in vacuum [28]. 
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